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ABSTRACT 

Observations of radio halos and relics in galaxy clusters indicate efficient elec- 
tron acceleration. Protons should likewise be accelerated and, on account of weak 
energy losses, can accumulate, suggesting that clusters may also be sources of very 
high-energy (VHE; E > 100 GeV) gamma-ray emission. We report here on VHE 
gamma-ray observations of the Coma galaxy cluster with the VERITAS array 
of imaging Cherenkov telescopes, with complementing Fermi-LAT observations 
at GeV energies. No significant gamma-ray emission from the Coma cluster was 
detected. Integral flux upper limits at the 99% confidence level were measured 
to be on the order of (2 - 5) x 10~ 8 ph. m' 2 s" 1 (VERITAS, > 220 GeV) and 
~2x 10 -6 ph. mr 2 s^ 1 (Fermi, 1 — 3 GeV), respectively. We use the gamma-ray 
upper limits to constrain CRs and magnetic fields in Coma. Using an analytical 
approach, the CR-to-thermal pressure ratio is constrained to be < 16% from 
VERITAS data and < 1.7% from Fermi data (averaged within the virial radius). 
These upper limits are starting to constrain the CR physics in self-consistent cos- 
mological cluster simulations and cap the maximum CR acceleration efficiency at 
structure formation shocks to be < 50%. Alternatively, this may argue for non- 
negligible CR transport processes such as CR streaming and diffusion into the 
outer cluster regions. Assuming that the radio-emitting electrons of the Coma 
halo result from hadronic CR interactions, the observations imply a lower limit 
on the central magnetic field in Coma of ~ (2 — 5.5) fiG, depending on the radial 
magnetic-field profile and on the gamma-ray spectral index. Since these values 
are below those inferred by Faraday rotation measurements in Coma (for most of 
the parameter space), this renders the hadronic model a very plausible explana- 
tion of the Coma radio halo. Finally, since galaxy clusters are dark-matter (DM) 
dominated, the VERITAS upper limits have been used to place constraints on 
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the thermally-averaged product of the total self-annihilation cross section and 
the relative velocity of the DM particles, iav) . 

Subject headings: cosmic rays — gamma rays: galaxies: clusters — galaxies: clusters: 
general — galaxies: clusters: individual: Coma (ACO 1656) — dark matter - 
magnetic fields 



-6 - 



Introduction 



Clusters of galaxies are the largest virialized objects in the Universe, with typical sizes 
of a few Mpc and masses on the order of 10 14 to 1O 15 M . According to the currently favored 
hierarchical model of cosmic structure formation, larger objects formed through successive 
mergers o f smaller objects with galaxy clusters sitting on top of this mass hierarchy (see 



Voit 



2005 



for a review). Most of the mass (~80%) in a cluster is dar 
indicated by galaxy dynamics and gravitational lensing (IDiaferio et al. 



t mat ter (DM), as 



20081 ). Baryonic gas 



making up the intra-cluster medium (ICM) contributes about 15% of the total cluster mass 
and individual galaxies account for the remainder (about 5%). The ICM gas mass also 
comprises a significant fraction of the observable (baryonic) matter in the Universe. 



The ICM is a hot (T ~ 10 8 
X-ray regime (see, e.g., iPetrosian 



'£) pla sma emitting thermal bremsstrahlung in the soft 



20011 ). This plasma has been heated primarily through 
collisionless structure-formation shocks that form as a result of the hierarchical merging 
and accretion processes. Such shocks and turbulence in the ICM gas in combination with 



intr a-cluster magnetic fields also 



e.g. 



Colafrancesco fc Blasi 



1998 



provide a mean s to accelerate particles efficiently (see, 



Ryu et al. 



20031 ). Many clusters feature megaparsec scale 



halos of nonthermal radio emission, indic ative of a population of relativistic electrons 



and magnetic fields permeating the ICM ( jCassano et al. 



20101 ) . There are two competing 



theories to explain radio halos. In the "hadronic model" , the radio-emitting electrons and 
positrons ar e produced in inelastic collisions o f cosmic-ray (CR) ions with the thermal gas 



of the ICM (IDennison 



1980 



Enfilin et al 



20111 ) . In the "re-acceleration model", a long-lived 



pool of 100-MeV electrons — previously accelerated by formation shocks, galactic winds, or 
jets of active galactic nuclei (AGN) — interacts with plasma waves that are excited during 
states of strong ICM turbulence, e.g., after a cluster merger. This may result in second 
order Fermi acceleration and may produce energetic electrons (~ 10 GeV) sufficient to 
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explain the observable radio emission (ISchlickeiser et al.lll987t iBrunetti fc Lazarianll201ll ) 



Observations o 



Sarazin fc Lieu 



2004 



possi bly nonthermal emi ssion from clusters in the extreme ultraviolet (EU V 



Eckert et al. 



1983) and hard X-rays (IRephaeli &: Gruber 



2002 



Fusco-Femiano et al. 



20071 ) may provide further indication of relativistic particle populations 



in clusters, although the interpretation of these observations as nonther mal diffuse emissio n 



has been disputed on the basis of more sensitive observations (see, e.g., lAjello et al. 

ft 



2009 



201 



Wiketal 



20091 ) 



Galaxy clusters have, for many years, been proposed as sources of gamma rays. If shock 
acceleration in the ICM is an efficient process, a population of highly relativistic CR protons 
and heavy ions is to be expected in the ICM. The main energy-loss mechanism for CR 
hadrons at high energies is pion production through the interaction of CRs with nuclei in the 
ICM. Pions are short lived and decay. The decay of neutral pions produces gamma rays and 
the decay of charged pions produces muons, which then decay to electrons and positrons. 
Due to the low density of the ICM (tiicm ~ 10~ 3 cm -3 ), the large size and the volume-filling 
magnetic fields in the ICM, CR hadrons will be confined in the cluster on timescales 



1996 



Berezinsky et al. 



19971 ) 



comparable to, or longer than, the Hubble time (IVolk et al. 
and they can therefore accumulate. For a given CR distribution function, the hadronically 
induced gamma-ray flux is dire ctly proportional t o the CR-to-thermal pressure fraction, 



^cr = (Pgr) I <Pth> (see, e.g. 



Enfilin et al. 



20071 ). where the brackets indicates volume 



avera ges. A very modest o f a few percent implies an observable flux of gamma rays 



(.e.g. 



Pfrommer fc Enfilin 



2004ah . 



Hydrostatic estimates of cluster masses, which are determined by balancing the thermal 
pressure force and the gravitational force, are biased low by the presence of any substantial 
nonthermal pressure component, including a CR pressure contribution. Similarly, a 
substantial CR pressure can bias the temperature decrement of the cosmic microwave 
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background (CMB) due to the Sunyaev-Zel'dovich effect in the direction of a galaxy cluster. 
This could then severely jeopardize the use of clusters to determine cosmological parameters. 
Comparing X-ray and optical potential profiles in the centers of galaxy clusters yields an 
upper limit of 20-30% of nonthermal pressure (that can be composed of CRs, magn etic 



2008 



20101 ). An 



fields or turbulence) relative to the thermal gas pressure (IChurazov et al.l 
analysis that compares spatially resolved weak gravitational lensing and hydrostatic X-ray 
masses for a sample of 18 galaxy clusters detects a deficit of the hydrostatic mass estimate 
compared to the lensing mass of 20% at -R500 - the radius within which the mean density is 
500 times the critical density of the Universe - suggesting again a substantial nonthermal 



pressure contribution on large scales (jMahdavi et al. 



20081 ) . Observing gamma-ray emission 



is a complementary method of constraining the pressure contribution of CRs that is most 
sensitive to the cluster core region. However, it assumes that the CR component is fully 
mixed with the ICM and may not allow for a detection of a two-phase structure of CRs 
and the thermal ICM. An Xqr of only a few percent is required in order to produce a 
gamma-ray flux observable with the current generation of gamma-ray telescopes, rendering 
this technique at least as sensitive as the dynamical and hydrostatic methods (which are 
more general in that they are sensitive to any nonthermal pressure component). 

Gamma-ray emission can also be produced by Compton up-scattering of ambient 
photons, for example CMB photons, on ultra-relativistic electrons. Those electrons can 
either be secondaries from the CR interactions mentioned above, or injected into the ICM 
by powerful cluster members and further accelerated by d iffusive shock acceleration or 



turbulent reacceleration processes ( ISchlickeiser et al. 



19871 and references therein). 



A third mechanism for gamma-ray production in a galaxy cluster could be self- 
annihilation of a DM particle, e.g., a weakly interacting massive particle (WIMP). As 
already mentioned, about 80% of the cluster mass is in the form of dark matter, which makes 
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galaxy clusters interesting targets for DM se arches (jEvans et al.ll2004l ; iBergstrom fc Hooper 



2006 



Pinzke et al 



2009 



Cuesta et al. 



201 if ) despite their large distances compared to 



other common targets fo r DM searches, s uch as dwarf spheroidal g alaxies (IStrigari et al. 



2007 




Acciari et al. 


20 


10; 


Alii 


i et al. 


2009 


) or the Gal 


Aharonian et al. 


2006, 


2009a: 


Abramowski et al. 


20111 



actic Center ( Kosack et al. 



2004J ; 



While several observations of clusters of galaxies have been made with satellite-borne 



and ground-based gamma-ray telescopes, a detec tion of gamma-ray em i ssion from a c 



has yet to be made. Observations with EGRET ( jSreekumar et al. 



1996 



Reimer et al. 



uster 



20031 ) 



and the Large Area Te 



(Ackermann et al. 



escope (LAT) on board the Fermi Gamma-ray Space Telescope 
2010t ) have provided upper limits on the gamma-ray fluxes (typically 
~ 10~ 9 ph cm 2 s _1 for Fermi-LAT observations) for several galaxy clusters in the MeV 
to GeV band. Upper limits on the very-high-energy (VHE) gamma-ray flux from a small 
sample of clusters, including the Coma cluster, have been provid ed by observations with 



ground-based imaging a 


tmospheric Cherenkov telescopes 


flACTs; 


Perkins et al. 


2006 


Aharonian et al. 


2009b: 


Aleksic et al.l 


201o|; 


Aleksic et al. 


2012 


)■ 





2008 



The Coma cluster of gala xies (ACQ 1656) i s one of the most thoroughly studied 



clusters acr oss all wavelengths (|Voges et al. 



(z = 0.023; 



Smith 



Struble & Rood 



1998 



Kubo et al 



Thierbach et al. 



19991 ). Located at a distance of about 100 Mpc 
20031 ) . it is one of the closest massive clusters (M ~ 10 15 M p ; 



20081 ). It hosts both a giant radio halo ( jGiovannini et al. 



1993 : 



20031 ) and peripheral radio reli c, which appears connec ted to the radio 



halo with a "diffuse" bridge (see discussion in 



Brown & Rudnick 



suggested (jEnfilin et al. 



20111 ). It has been 



19981 ) and successively d emonstrated by cosmological simu 



which model the non t hermal emission pro cesses (IPfrommer et al. 



2008b 



Pfrommer 



ations 



2008 



Battaglia et al. 



2009 



Skillman et al 



20111 ). that the relic could well be an infall shock. 



Extended soft thermal X-ray (SXR) emission is evident from the ROSAT all-sky survey in 
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the .1 to 2.4 keV band flBriel et al.lll992f ). Observations with XMM- Newton f lBriel et al. 



20011 ) revealed substructure in the X-ray halo supported 



of at least ~ 10% of the total pressure ( jSchuecker et al. 



jy su bstantial turbulent pressure 



20041 ). The Coma cluster is a 



natural candidate for gamma-ray observations. 



In this article, results from the VERITAS observations of the Coma cluster of galaxies 
are reported, with complementing analysis of available data from the Large Area Telescope 
(LAT) on board the Fermi Gamma-ray Space Telescope. The VERITAS and Fermi-LAT 
data have been used to place constraints on cosmic-ray particle populations, magnetic fields, 
and dark matter in the cluster. Throughout the analyses, a present day Hubble constant of 
H = 100/i km s _1 Mpc -1 with h = 0.7 has been used. 



2. VERITAS Observations, Analysis, and Results 



The VERITAS gamma-ray detector (1 Weekes et al 



2002 ) is an array of four 12 



Holder et al. 



20061 ) located at an 



m-diameter imaging atmospheric Cherenkov telescopes (1 
altitude of ~1250 m a.s.l. at the Fred Lawrence Whipple Observatory in southern Arizona 
(31° 40' 30" N, 110° 57' 07" W). Each of the telescopes is equipped with a 499-pixel camera 
covering a 3.5° field of view. The array, completed in the fall of 2007, is designed to detect 
gamma-ray emission from astrophysical objects in the energy range from 100 GeV to more 
than 30 TeV. Depending on the zenith angle and quality selection criteria imposed during 
the data analysis, the effective energy range may be narrower than that. The energy 
resolution is ~ 15% and the angular resolution (68% containment) is ~ 0.1° per event at 
1 TeV and slightly larger at low energy. At the time of the Coma cluster observations, 
the sensitivity of the array allowed for detection of a point source with a flux of 1% of the 
steady Crab Nebula flux above 300 GeV at the confidence level of five standard deviations 



(5a) in under 45 hours 

The Coma cluster was observed with VERITAS between March and May in 2008 
with all four telescopes fully operational. The total exposure amounts to 18.6 hours 
of quality-selected live time, i.e., time periods of astronomical darkness with clear sky 
conditions and no technical problems with the array. The center of the cluster was tracked 
in wobble mode, where the expected source location is offset from the center of the field o f 



view by 0.5 degrees, to allow for simultaneous background estimation (jFomin et al. 



19941 ) 



All of the observations were made in a small range with average zenith angle ~ 21°. 



in 



The data analysis was performed fo 



Cogan et al. 



(120071) and 



Daniel et al- 



l owing the standard VERITAS procedures described 



( 120071 ). Prior to event reconstruction and selection, 



all shower images are calibrated and cleaned. Showers are then reconstructed for events 
with at least two telescopes contributing images that pass the following quality selection 
criteria: more than four participating pixels in the camera, number of photoelectrons in the 
image larger than 75, and the distance from the image centroid to the center of the camera 
less than 1.43°. These quality selection criteria impose an energy thresholds of about 220 
GeV. In addition, events for which only images from the two closest-spaced telescopes^ 
survive quality selection are rejected, as they introduce an irreducible high background rate 



x The integral flux sensitivity above 300 GeV was improved by about 30% with the relo- 
cation of one telescope in the summer of 2009. 

2 The energy threshold is defined as the energy corresponding to the maximum of the 
product function of the observed spectrum and the collection area. It does not vary sig- 
nificantly for the different source scenarios and assumed spectral indices reported in this 
work. 

3 In the array configuration prior to summer 2009, two telescopes had a separation of only 
35 m. 
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due to local muons, degrading the instrument sensitivity (IMaier fc Knappl 120071 ). 



Gamma-ray-like events are separated from the CR backgrou nd by imposing selectio n 



criteria (cuts) on 



Krawczynski et al. 
images iHillasl 



re m ean-scaled length and width parameters (lAharonian et al. 



1997 



20061 ) calculated from a parametrized moment analysis of the shower 
19851 ). These parameters are averages over the four telescopes weighted with 
the total amplitude of the images, that measure the image moments width and length 
scaled with values expected for gamma rays. In this analysis, events with a mean-scaled 
length in the range 0.05-1.19 and a mean-scaled width in the range 0.05-1.08 are selected as 
gamma-ray-like events. These ranges for the gamma-hadron separation cuts were optimized 
a priori for a weak point source (3% Crab Nebula flux level) and a differential spectral index 
of 2.4, using data taken on the Crab Nebula during the same epoch. Because the VHE 
gamma-ray spectru m for the Coma cluster is expected to be a power-law function with an 



index of about 2.3 ( IPinzke &: Pfrommer 



20101 ). these cuts are suitable for the analysis of 



the Coma cluster data set. It is noted that slightly varying the spectral index (± 0.2) does 
not significantly impact the cuts used for quality selection and gamma-hadron separation 
in this work. 

The Coma cluster is a very rich cluster of galaxies with many plausible sites for 
gamma-ray emission: the core region, the peripheral radio relic, and individual powerful 
cluster member galaxies. VERITAS has a large enough field of view to allow investigation 
of several of these scenarios. In this work, the focus has been on the core region and 
three cluster members. The core region is treated as either a point source or a mildly 
extended source, a uniform disk with intrinsic radius 0.2° or 0.4°, similar to the extension 
of the thermal soft X-ray emission from the core. There is eviden ce of a recent merger 



event between the two central galaxies NGC 4889 and NGC 4874 ( ITribble 



19931 ) . There 



is also evidence for an excess of nonthermal X-ray emission from these galaxies as well as 



13 



from the galaxy NGC 4921 (INeumann et al.ll2003l ). Therefore, searches for point-like VHE 
gamma-ray emission have been conducted at the locations of these galaxies. The regions of 
interest considered in this work are summarized in Table [TJ 



The ring-background model (lAharonian et al.ll200ll ) is used to estimate the background 
due to CRs misinterpreted as gamma rays (the cuts described above reject more than 99% 
of all CRs). The total number of events in a given region of interest is then compared to 
the estimated background from the off-source region scaled by the ratio of the solid angles 
to produce a final excess o r deficit. The V HE gamma-ray significance is then calculated 



according to Formula 17 in 



Li &: Mai (Il983l ). Significance skymaps over the VERITAS field 



of view produced with a 0.2° integration radius are s 



rown in Figure CD with overlaid X-ray 



and radio co ntours from the RO SAT all-sky survey (jBriel et al. 



observations ( jBrown &: Rudnick 



19921 ) and GBT 1.4 GHz 



201 ll ) respectively. 



Depending on the assumed extent of the source and the point-spread function, we can 
define an ON region, into which a defined fraction of the source photons should fall. No 
significant excess of VHE gamma rays from the Coma cluster was detected with VERITAS, 
as illustrated by the 9 2 distribution shown in Figure [2j in which source events would pile up 
at small values of 9 2 for a point source and fall into a somewhat wider range of 9 2 values for 
an extended source. The 9 2 distribution is a plot of event density versus the square of the 
angular separation from a given location. It permits a comparison of the ON-source event 
distribution with that of other locations, in this case a ring-shaped region, into which only 
background events should fall, the so-called OFF-source region. The 9 2 distribution extends 
out to 0.42 square degrees to cover both the case of point-like and extended emission from 
the core of the Coma cluster. The 9 2 distributions for the member galaxies also considered 
in this work are very similar to that in Figure [2] and show no excess of gamma rays. A 
99% confidence level upper limit is calculated for each region of interest using events from 
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the ON-source and OFF-source regions and the method described by iRolke et al.1 ( 120051 ) 
assuming a Gaussian-distributed background. A lower bound of zero is imposed on the 
gamma-ray flux from the Coma cluster, which prevents artificially low flux upper limits in 
the case that the best-fit source flux is formally negative. Figure [3] shows the distribution of 
significances over the VERITAS skymap, which is well fit by a Gaussian with a mean close 
to zero and a standard deviation within a few percent of unity. 

Table [2] lists the upper limits for the selected regions of interest shown in Table [TJ 
These upper-limit calculations depend on the gamma-ray spectrum, which in this work 
is assumed to be a power law in energy, dN/dE oc E~ a , where the spectral index a was 
allowed to have a value of 2.1, 2.3, or 2.5. 



3. Fermi-LAT Analysis and Results 

LAT on board Fermi has observed the Coma cluster in all-sky survey mode since 
its launch in June 2008. Fermi-LAT is sensitive to gamma rays in the 2 MeV to ~ 300 
GeV energy range and is complementary to the VERITAS observations. 



Ackermann et al. 



( 120101 ) reported on the search for gamma-ray emission from thirty-three galaxy clusters in 
the data from the first 18 months, including the Coma cluster, for which an upper limit 
of 4.58 x 10~ 9 ph cm -2 s -1 in the 0.2 to 100 GeV energy band was reported. This limit 
is expected to improve as the exposure is increased. In this work an updated analysis is 
presented as a complement to the VERITAS results which includes data taken between 
August 5, 2008 and April 17, 2012. 

The LAT- d ata a nalysis of this work follows the same procedure as described in detail 



111 



Nolan et al. 



(120121 ) and was performed with the Fermi Science Tools version 9.23.1. To 



only include events with high probability of being photons, the P7SOURCE class and the 
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corresponding P7SOURCE_V6 instrument-response functions were used throughout this 
work. 

A zenith-angle cut of 100° was applied to eliminate albedo gamma rays from the 
Earth's limb, excluding time intervals during which any part of the region of interest (ROI) 
was outside the field of view. In addition, time intervals were removed during which the 
observatory was transiting the Southern Atlantic Anomaly or the rocking angle exceeded 
52°. 

The ROI is defined to be a square region of the sky measuring 14° on a side and 
centered on aj 2 ooo — 194.953 and £j2ooo — 27.9806, the nominal position of the Coma 
cluster. 

Only photons with reconstructed energy greater than 1 GeV are considered, for which 
the 68%-containment radius of the point-spread function (PSF) is narrower than ~ 0.8°. 
The Fermi-LAT collaboration estimates the systematic uncertainties on the effective area 
at 10 GeV to be around 10% Q. 



The background emission in t he ROI was mode 



in the second LAT source catalog (jNolan et al 



ed using fourteen point sources listed 



20121 ). the LAT standard Galactic diffuse 



emission component (gal_2yearp7v6_v0 . f it), and the corresponding isotropic template 
(iso_p7v6source.txt) that accounts for extragalactic emission and residual cosmic-ray 
contamination. Due to the large tails of the PSF at low energy, further fourteen point 
sources, lying ~ 4° outside the ROI, were included in the source model. 

The energy spectra of twenty-four sources are described by a power law. The remaining 
four sources.], being bright sources, are modeled with additional degrees of freedom using 



4 http: / /fermi. gsfc.nasa.gov/ssc/data/analysis/LAT_caveats. html 

5 2FGLJ1303. 1+2435, 2FGLJ1310.6+3222, 2FGLJ1226.0+2953, and 2FGLJ1224.9+2122 
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the log-normal representation, which is typically used for modeling Blazar spectra. 



The analysis is performed in three energy bins: 1-3 GeV, 3-10 GeV, and 10-30 GeV. To 



find t he best fit spectral parameters, a binned maximum-likelihood analysis ( jMattox et al. 



19961 ) is performed for each energy bin on a map with 0.1° pixel size in gnomonic (TAN) 
projection, covering the entire ROI. To determine the significance of the sources, and in 
particular that of th e Coma cluster, the analysis tool uses the likelihood-ratio test statistic 



(jMattox et al. 



19961 ) defined as, 



TS = -2(lnL -InL), (1) 

where L is the maximum likelihood value for the null hypothesis and L is the maximum 
likelihood with the additional source at a given position on the sky. 

In the likelihood analysis the spatial parameters of the sources were kept fixed at the 
values given in the catalog, whereas the spectral parameters of the point sources in the 
ROI, along with the normalization of the diffuse components, were allowed to freely vary. 
We analyzed three cases in which the gamma-ray emission from the Coma cluster was 
assumed to follow a power-law spectrum with a photon index a = 2.1, 2.3, and 2.5. The 
spectral indices of all point sources were permitted to freely vary between a = and a = 5. 
We considered the emission as being caused both by a point-like and a spatially extended 
source (a uniform disk) with radius r = 0.2° or r = 0.4°, as in the VERITAS analysis. 

No significant gamma-ray signal was detected. For one free parameter, the flux from 
the Coma cluster, the detection significance is computed as the square root of the test 
statistic (TS follows a \\ distribution). The highest test statistic was obtained for the 
high-energy band, where TS ~ 0.8 for the point source model, TS ~ 0.7 for the disk model 
with r = 0.2°, and TS ~ 2 for the disk model with r = 0.4°. 



We therefore used the profile likelihood method ( iRolke et al. 



20051) to derive flux 
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upper limits at the 99% confidence level in the energy range 1-30 GeV, assuming both an 
unresolved, point-like or spatially extended emission, as shown in Table [3J 



4. Gamma Ray Emission from Cosmic Rays 

We decided to adopt a multifaceted approach to constrain the CR-to-thermal pressure 
distribution in the Coma cluster using the upper limits derived from the VERITAS and 
Fermi-LAT data in this work. This approach includes (1) a simplified multi-frequency 
analytical model that assumes a constant CR-to-thermal energy density and a power-law 
spectrum in momentum, (2) an analytic model derived from cosmological hydrodynamical 
simulations of the formation of galaxy clusters, and (3) a model that uses the observed 
intensity profile of the giant radio halo in Coma to place a lower limit on the expected 
gamma-ray flux in the hadronic model - where the radio-emitting electrons are secondaries 
from CR interactions and which is independent of the magnetic field distribution. This last 
approach translates into a minimum CR pressure which, if challenged by tight gamma-ray 
limits/detections, permits scrutiny of the hadronic interaction model of the formation of 
giant radio halos. Alternatively, realizing a spatial CR distribution that is consistent with 
the flux upper limits, and requiring the model to match the observed radio data, enables us 
to derive a lower limit on the magnetic field distribution. We stress again that this approach 
assumes the validity of the hadronic interaction model. Modeling the CR distribution 
through different techniques enables us to bracket our lack of understanding about the 
underlying plasma physics that shapes the C R distribution hen ce to reflect the Bayesian 



priors that are imposed on the modeling (see 



Pinzke et al. 



2011 



, for a discussion). 
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4.1. Simplified analytical model 



We start by adopting a simplified analytical model that assumes a power-law CR 
spectrum and a constant CR-to-t hermal pressure ratio , i.e., w e adopt the isobaric model 



of CRs following the approach of 



Pfrommer &: Enfilinl (j2004al ). To be independent of 



additional assumptions and in line with earlier work in the literature, we do not impose 
a low-momentum cutoff, q, on the CR distribution function, i.e., we adopt q = 0. Since, 
a priori, the CR spectral index is unconstrained^ we vary it in the range 2.1 < a < 2.5, 
which is compatible with the radio spectral index of the giant radio halo of the Coma 



cluster after accounting fo r the spectra 



Sunyaev-Zel'dovich effect (jEnfilin 



2002 



steepening at frequencies v ~ 5 GHz due to the 



Pfrommer fc Enfilin 



2004af )J I To model the thermal 



pressure, we adopt the electron de nsity profile for t he Coma cluster that has been inferred 



from ROSAT X-ray observations (IBriel et al. 



1992J) and use a constant temperature of 



kT = 8.25 keV throughout the virial region. 

Table H] shows the resulting constraints on the CR-to-thermal pressure ratio, 
A^cr = (Pcr) I (Pth), averaged within the virial radius, i? vir = 2.2 Mpc, that we define as 
the radius of a sphere enclosing a mean density that is 200 times the critical density of 
the Universe. Constraints on Xqr with VERITAS flux upper limits (99% CL) strongly 
depend on a. This is due to the comparably large energy range from GeV energies (that 
dominate the CR pressure, provided a > 2 and the CR population has a nonrelativistic 



6 The hadronic interaction physics guarantees that the CR spectral index coincides with 
that of the resulting pion-decay gamm a-ray emission at energies E 3> 1 GeV that are well 



above the pion bump (see discussion in 



Pfrommer fc EnBlin 



2004aJ) 



7 Assuming a magnetic field of 1 /iG, the CR protons responsible for the GHz radio 
emitting electrons have an energy of ~ 100 GeV and are ~ 20 times less energetic than those 
CR protons responsible for 200-GeV gamma-ray emission. 
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low- momentum cutoff, i.e., q < m p c, where m p is the proton mass) to energies at 220 
GeV, where our quality selection criteria imposed the energy threshold. These gamma-ray 
energies correspond to 1.6 TeV CRs - an energy ratio of more than 3 orders of magnitude, 
which explains the sensitivity to small changes in a. The flux measurements within 0.2° 
are the most constraining due to a competition between the integrated signal and the 
background as the integration radius increases. This yields limits on Xcr between 0.048 
and 0.43 (for a varying between 2.1 and 2.5), with a constr aint of Xrn < 0.1 for a = 2 .3 



Pinzke fc Pfrommeri (120 101 ) 



(close to the spectral index predicted by the simulations of 
around 220 GeV). Constraints on Xcr with Fermi-LAT limits (99% CL) depend only 
weakly on a because GeV-band gamma rays are produced by CRs with energies near the 
relativistic transition, that dominantly contribute to the CR pressure. XcR-constraints with 
Fermi-LAT limits are most constraining for an aperture of 0.4°; despite the slightly weaker 
flux upper limits in comparison to the smaller radii of integration, we expect a considerably 
larger gamma-ray luminosity due to the increasing volume in this model. The best limit 
of Xqr < 0.012 is achieved for a = 2.3, while the limit for a = 2.1 is only slightly worse 
(X CR < 0.017). 



4.2. Simulation-based approach 



We complement the simplified analytical analysis with a more realistic and predictive 
approach derived from cosmological hydrodynamical s imulations. We adop t the u niversal 



Pinzke fc Pfrommeri (120101 ) to 



spectral and spatial gamma-ray model developed by 
estimate the emission from decaying neutral pions which in clusters dominates over 
the inverse-Compton (IC) emission above 100 MeV. Given a density profile as, e.g., 
inferred by cosmological simulations or X-ray observations, the analytic approach models 
the CR distribution and the associated radiative emission processes from radio to the 
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gamma-ray band. This formalism was derived from high-resolution simulations of 
clusters of galaxies that included radiative hydrodynamics, star formation and supernova 
feedback, and it followed the CR physics by tracing the most important injection and loss 
processes self-consistently while accounting for the CR pressure in the equation of motion 



( iPfrommer et al 



2006 



Enfilin et al 



2007 



Jubelgas et al 



20081 ). The results are in line 



with earlier numerical results on some of the overall characteristics of the CR distribution 



and the associated radiative emission proc e sses (IDolag fc EnBlin 



Miniati 



2003 



Pfrommer et al. 


2007, 


2008b: 


Pfrommer 


2008) 



2000 : 



Miniati et al 



2001 



The overall normalization of the CR and gamma-ray distribution scales nonlinearly 
with the acceleration eff iciency at structure formation shocks. Following re cent observatio ns 



of su pernova remnants (IHelder et al. 



20091 ) as well as theoretical studies (IKang fc Jones 
20051 ). we adopt an optimistic but nevertheless realistic value of this parameter and assume 
that 50% of the dissipated energy at strong shocks is injected into CRs, with this efficiency 
decreasing rapidly for weaker shocks. Since the vast majority of internal for mation shock s 



Ryu et al. 



(merg er and flow shocks) are weak shocks with Mach numbers M < 3 (e.g. 
20031 ). they do not contribute significantly to the CR population in clusters. Instead, 
strong shocks during the formation epoch of clusters and strong accretion shocks at the 
present time (at the boundary of voids and filaments/supercluster regions) dominate the 
acceleration of CRs which are adiabatically transported through the cluster. Hence, the 
model provides a plausible upper limit for the CR contribution from structure formation 
shocks in galaxy clusters which can be scaled with the effective acceleration efficiency. 
Other possible CR sources, such as AGN and starburst- driven galactic winds have been 
neglected for simplicity but could in principle increase the expected gamma-ray yield. 

These cosmological simulations only consider advective transport of CRs by bulk 
gas flows that inject a turbulent cascade, leading to centrally-enhanced density profiles. 
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However, other means of CR transport such as diffusion and streaming may flatten the CR 
radial profiles. The CRs stream along magnetic field lines in the opposite direction of the 
CR number density gradient (at any energy). In the stratified cluster atmosphere, this 
implies a net flux of CRs towards larger radii, equalizing the CR number density with time 
if not counteracted by advective transport. It has been suggested that advection velocities 
only dominate over the CR streaming velocities for periods with trans- and supersonic 
cluster turbulence during a cluster merger and drop below the CR streaming velocities for 
relaxing clusters. As a consequence, a bimodality of the CR spatial distribution is expected 
to result; with merging; (relaxed) clusters showing a centrally concentrated (flat) CR energy 



density profile (jEnBlin et al. 



20111 ) . This translates into a bimodality of the expected diffuse 



radio and gamma-ray emission of clusters, since more centrally concentrated CRs will find 
higher target densities for hadronic CR proton interactions. As a result of th is, relaxed 



clust ers could have a reduced gamma-ray luminosity by up to a factor of five (jEnBlin et al. 



201ll ). Hence, tight upper limits on the gamma-ray emission can constrain a combination 
of acceleration physics and transport properties of CRs. 

We adopt the density profile of thermal electrons as discussed in §4.11 and model 
the temperature profile of the Coma cluster with a constant central temperature of 
kT = 8.25 keV and a characteristic decline toward the cluster periphery in accordance 
with a fit to the universa l profile derived from cosmological cluster simulations 



( jPinzke &: Pfrommer 



2010 



deep Chandra cluster data (IVikhlinin et al. 



'frommer et al. 



2007]) and the behavior of a nearby sample of 



20051 ). This enables us to adopt 



spectral distribution of CRs according to the model by 



; he spa tial and 



Pinzke & Pfrommer 



(|2010| ) that 



neglects the contribution of supernova remnants, AGN, and cluster galaxies. 

Figure H] shows the expected integral spectral energy distribution of Coma within 
the virial radius (dotted line). This suggests a spectral index of a = 2.1 in the energy 
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interval 1-3 GeV and a = 2.3 for energies probed by VERITAS (> 220 GeV). Also shown 
are integrals of the differential spectrum for finite energy intervals across the angular 
apertures tested in this study (dashed lines). These model fluxes (summarized in Table 
H]) are compared to Fermi and VERITAS flux upper limits for the same energy intervals. 
Constraints on Xcr with the gamma-ray flux limit of Fermi in the energy interval 1-3 GeV 
(< 0.4°) are most constraining, since that combination of a specific energy interval and 
aperture minimizes the ratio of the upper limit to the expected model flux. In particular, 
this upper limit is 24% below the model predictions that assume an optimistically large 
shock-acceleration efficiency and CR transport parameters as laid out above. Hence this 
enables us to constrain a combination of maximum shock acceleration efficiency and CR 
transport parameters. In our further analysis, we use the most constraining Fermi-LAT 
flux limits in the energy interval 1-3 GeV as well as the gamma-ray flux limits of VERITAS 
in the energy range above 220 GeV. 

Figure [5] shows the CR-to-thermal pressure ratio, X C r = (Pcr) / (Pth), as a function 
of radial distance, R, from the Coma cluster center and contained within R. All radii are 
shown in units of the virial radius, i? v i r = 2.2 Mpc. To compute the CR pressure, we assume 
a low-momentum cutoff of the CR distribution at q — 0.8 m p c, where m p is the proton mass. 
This is suggested by cosmological cluster simulations and reflects the high Coulomb cooling 
rates at low CR energies. The CR-to-thermal pressure ratio rises toward the outer regions 
on account of the higher efficiency of CR acceleration at the peripheral accretion shocks 
compared to the weak central flow shocks. Adiabatic compression of a mixture of CRs and 
thermal gas disfavors the CR pressure relative to the thermal pressure on account of the 
softer equation of state of CRs. The weak increase of Acr toward the core is due to the 
comparably fast thermal cooling of gas. 

In the case of VERITAS, for the most constraining regions tested (within an aperture 
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of radius 0.2°), the predicted CR pressure is a factor of 7.2 below the inferred upper 
limits of VERITAS (see Table |5] and assuming a spectral index of a = 2.3 which matches 
the simulated one at energies E 1 = 200 GeV). To first order, we can scale the averaged 
CR-to-thermal pressure ratio of our model by that factor, keep the spatial behavior and 
obtain an integrated limit of the CR-to-thermal pressure ratio of Xcr < 0.112 within 0.2° 
that translates to a limit within the cluster virial radius of Xcr < 0.162 (solid lines of Figure 
[5]). This limit is less constraining by 50% in comparison to the simplified analytical model, 
which gives Xcr < 0.1. This difference is explained by the concavity of the simulated 
spectrum which therefore carries more pressure at GeV energies than a pure power-law 
spectrum with a = 2.3. 

As already aluded to, the most constraining Fermi-LAT upper limit in the energy 
interval 1-3 GeV (< 0.4°) is a factor of 0.76 smaller than our model predictions (assuming 
a = 2.1 which is very close to the simulated spectral index for the energy range 1-3 GeV). 
Scaling our integrated CR-to-thermal pressure profile yields a constraint of Xcr < 0.012 
within 0.4° that translates to a limit within the cluster virial radius of Xcr < 0.017 
(dashed lines of Figure |5]). The X C r constraint evaluated within the cluster virial radius is 
comparable to the constraint of Xcr < 0.017 in our simplified model. Naturally, with the 
Fermi-LAT limits we probe the region around GeV energies that dominate the CR pressure, 
and we do not expect any differences to the simplified power-law model in comparison to 
the universal CR spectrum with its concave CR spectrum found in the simulations. 



4.3. Minimum gamma-ray flux 

For clusters that host radio halos, we can derive a minimum gamma-ray flux in the 
hadronic model of radio halos - where the radio-emitting electrons are secondaries from CR 
interactions. Hadronic interactions channel about the same power into secondary electrons 
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and 7r°-decay gamma rays. A stationary distribution of CR electrons loses all its energy to 
synchrotron radiation for strong magnetic fields (B S> -Bcmb — (1 + z) 2 3.2/iG, where -Bcmb 
is the equivalent magnetic field strength of the CMB so that B^, mb /8tc equals the CMB 
energy density). Thus the ratio of gamma-ray to synchrotron flux becomes independent of 



1989 



Pohl 



1994 



Pfrommer 



2008|), in 



the spatial distribution of CRs and thermal gas ( IVolkl 
particular with a v ~ 1 as the observed synchrotron spectral index. Hence we can derive a 
minimum gamma-ray flux in the hadronic model 



.4. 



7, mm 



A v AttD 



2 ' 
lum 



(2) 



where L v is the observed luminosity of the radio mini-halo, Di nm denotes the luminosity 
distance to the respective cluster, and A~ and A v are dimensional constants tha t depend on 



the hadronic physics of the interaction (jPfrommer 



2008 



Pfrommer et al 



2008aj). Lowering 



the magnetic field would require an increase in the energy density of CR electrons to 
reproduce the observed synchrotron luminosity and thus increase the associated gamma-ray 
flux. 

To derive a minimum gamma-ray flux that can be compared to the upper limits, we 
need to determine the radio flux within the corresponding angular regions. To this end, 
we fit the point-s ource-subtracted, azimuthally-averaged radio-halo profile at 1.38 GHz 



(IDeiss et al. 



19971 ) with a /3-model, 

S v (r±) = S 

where Sq = 1.1 x 10~ 3 Jyarcmin -2 , r c = 



Tc. 



-3/3+1/2 



(3) 



450 kpc, and = 0.78 . With in the error bars, this 



Govoni et al. 



(120011 ) when scaled with a 



profile is consistent with 326-MHz data taken by 
radio spectral index of 1.15. 

The results for the minimum gamma-ray flux -F 7iIIim (> 220 GeV) and the minimum 
CR-to-thermal pressure ratio X CRmin = X CR F 7jmin /F 7jiso are shown in Tabled where 
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F 7i i SO is the gamma-ray flux in the simplified model introduced in §4.11 Even in the most 



constraining cases, and assuming a < 2.3, these are a factor of ~ 60 below the VERITAS 
upper limits (for a = 2.1, < 0.2°) and a factor of ~ 20 below the Fermi-LAT upper 
limits (for a = 2.3, < 0.4°). Note that these minimum gamma-ray fluxes are sensitive 
to the variation of the CR proton spectral index with energy as a result of, for example, 
momentum-depen dent diffusion. Assum ing a plausible value for the central magnetic field 



of Coma of 5 fiG ( IBonafede et al 



20101 ). the radio-halo emission at GHz frequencies is 
dominated by electrons with energy E e ~ 2.5 GeV (which corresponds to proton energies 
E p ~ 40 GeV). Gamma rays with an energy of 200 GeV are produced by CR protons 
with an energy of E p ~ 1.6 TeV - a factor of 50 higher than those probed by radio halo 
observations. A steepening of the CR proton spectral index of 0.2 between 40 GeV and 1.6 



TeV would imply a decrease in the minimum gamma-ray flux by a factor of two. 



4.4. Constraining the Magnetic Field 

In the previous section, we have obtained an absolute lower limit on the gamma-ray 
emission in the hadronic model by assuming high magnetic fields, B ^> -Bcmb- We can 
turn the argument around and use our upper limit on the gamma-ray emission (and by 
extension on the CR pressure) to infer a lower limit on the magnetic field needed to explain 
the observed radio emission. This, again, assumes the validity of the hadronic model of 
radio halos, in which the radio-emitting electrons are secondaries from CR interactions. A 
stronger gamma-ray constraint will tighten the magnetic-field limit. In case of a conflict 
with magnetic field measurements by other methods, e.g., Faraday rotation measure 
the hadronic model of radio halos would be challenged. The method we use to constrain the 

8 Generally, Faraday RM analyses of the magnetic field strength by, e.g., background 
sources observed through clusters, are degenerate with the magnetic coherence scale and 
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magnetic field inherits a dependence on the assumed radial scaling which we parametrize as 

'n e (r) 



Br) 



Bn 



n e (0) 



(4) 



as suggested by Faraday RM studies and numerical m agnetohydrodynamical (MHD) 



simulations (IBonafede et al. 



2010 



Bonafede et al 



denotes the Coma electron density profile (IBriel et al. 



2011 



and references therein). Here n e 



19921 ). In fact, the magnetic field 



in the Coma cluster is among the best constrained, because its proximity permits RM 
observations of seven radio sources located at projected distances of 50 to 15 00 kpc from the 



cluster center. The best-fit model yields B = 4.7^o; 8 /zG and olb = O.Slga (IBonafede et al 
2010l ). We aim to constrain the central field strength, Bo, and we permit the magnetic 
decline, ag, to vary within a reasonable range of = 0.2 as suggested by those Faraday 
RM studies. We proceed as follows: 



1. Given a model for the magnetic field with as and an initial guess for Bo, we 
determine the profile of the CR-to-thermal pressure ratio, XcR.(r), by matching the 
hadronically-produced synchrotron emission to the observed radio-halo emission over 
the entire extent. To this end, we deproject the fit to the surface-brightness profile o f 



Eq. [3] (using an Abel integral equation, see Appendix of 
yielding the radio emissivity, 



Pfrommer fc Enfilin 



2004aJ) 



Mr) 



6/3 



B 



3/3 = j„o (1 + r 2 /r c 2 ) 



-3/3 



(5) 



2vrr c (l+r 2 /r c 2 ) 3/3 ~ V2 
where B denotes the beta function. It is generically true for weak magnetic fields 
(B < -Bomb) m the outer parts of the Coma halo that the product X CR (r)X B (r) 
(where Xb denotes the magnetic-to-thermal energy density ratio) has to increase by a 



may be biased by the unknown correlation between magnetic and density fluctuations. 
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factor of about 100 toward the radio-halo periphery to account for the observed extent. 
If we were to adopt a steeper magnetic decline such as«B = 0.5 which produces a flat 
Xb(t), the CR-to-thermal pressure ratio would have to rise accordingly by a factor of 



2. Given this realization for Xcr, we compute the pion-decay gamma-ray surface- 
brightness profile, integrate the flux within a radius of (0.13,0.2,0.4) degree, and 
scale the CR profile in order to match the corresponding VERITAS/Ferrm flux upper 
limits. This scaling factor, X CR0 , depends on the CR spectral index, a, (assuming a 
power-law CR population for simplicity), the radial decline of the magnetic field, ag, 
and our initial guess for B . 

3. We then solve for B while matching the observed synchrotron profile and fixing 
the profile of Xcb.{t) as determined through the previous two steps. Note that for 
B 3> -Bomb and a radio spectral index of a v = 1, the solution would be degenerate 
since the luminosity of the radio halo scales as 



where Q(E) denotes the electron source function. 

4. Inverse- Compton cooling of CR electrons on CMB photons introduces a characteristic 
scale of -Bcmb — 3.2 /iG which imprints as a nonlinearity on the synchrotron emissivity 
as a function of magnetic field strength (see Eq. (jSJ)). Hence we have to iterate 
through the previous steps until our solution for the minimum magnetic field B Q 
converges. 



100. 





Table |6] shows the resulting lower limit of the central magnetic field ranging from 
B = 0.5 to 1.4 fiG in case the of VERITAS and from B Q = 1.4 to 5.5 /iG in the case 
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of Fermz-LATfj Since these lower limits on Bq are below the values favored by Faraday 
RM for most of the parameter space spanned by olb and a (and never exceed the values 
for the phenomenological Faraday RM- inferred 5-model), the hadronic model is a viable 
explanation of the Coma radio halo. In fact, the Fermi- LAT upper limits start to rule out 
the parameter combination of «b > 0.7 and a > 2.5 for the hadronic model of the Coma 
radio halo. Future gamma-ray observations of the Coma cluster may put more stringent 
constraints on the parameters of the hadronic model. 

A few remarks are in order. (1) For the VERITAS limits, the hardest CR spectral 
indices correspond to the tightest limits on B , because the CR flux is constrained around 
1 TeV and a comparably small fraction of CRs at 100 GeV would be available to produce 
radio-emitting electrons. A high magnetic field would be required to match the observed 
synchrotron emission. The opposite is true for the Fermi upper limits at 1 GeV, which 
probe CRs around a pivot point of 8 GeV: a soft CR spectral index implies a comparably 
small fraction of CRs at 100 GeV and hence a strong magnetic field is needed to match 
the observed synchrotron flux. (2) For a steeper magnetic decline (larger o>b), the CR 
number density needs to be larger to match the observed radio-emission profiles, which 
would yield a higher gamma-ray flux so that the upper limits are more constraining. This 
implies tighter lower limits for Bq. (3) Interestingly, in all cases, the 0.4°-aperture limits 
are the most constraining. For a given magnetic realization, a substantially increasing 
CR-to-thermal pressure profile is needed to match the observed radio profiles, and therefore 
that CR realization produces a larger flux within 0.4° in comparison with the simplified CR 
model (Xcr = const.), for which the 0.2°-aperture limits are more constraining in the case 
of VERITAS. Physically, the large CR pressure in the cluster periphery may arise from CR 



9 Note that a central magnetic field of 3 fiG corresponds in the Coma cluster to a magnetic- 
to-thermal energy density ratio of Xb = 0.005. 
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streaming into the large available phase space in the outer regions. 

As a final word, in Table M we show the corresponding values for the CR-to-thermal 
pressure ratio (at the largest emission radius at 1 Mpc) such that the model reproduces 
the observed radio surface-brightness profileo They should be interpreted as upper limits 
since they are derived from flux upper limits. For the Fermi-LAT upper limits, they range 
from 0.08 to 0.27; hence the Aqr profiles always obey the energy condition, i.e., Pgr < Pth, 
over the entire range of the radio-halo emission (< 1 Mpc)0 The corresponding values for 
Xqr in the cluster center are smaller than 0.01 for the entire parameter space probed in 
this study. We conclude that the hadronic model is not challenged by current Faraday RM 
data and is a perfectly viable possibility in explaining the Coma radio-halo emission. 



Emission from Dark Matter Annihilations 



As already mentioned in the introduction, most of the mass in a galaxy cluster is 
in the form of DM. While the nature of DM remains unknown, a compelling theoretical 



10 Note that in this section, we determine the radial behavior of Acr by adopting a specific 
model for the magnetic field and requiring the modeled synchrotron surface-brightness profile 
to match the observed data of the Coma radio halo. This is in contrast to the simplified 
analytical CR model where Acr is constant (§ 14. li and to the simulation-based model where 
Acr(?") is derived fr om cosmological clus t er simu lations (§ 14. 2p . 



ii 



See Figure 3 in 



minimum energy conditions, and 



Pfrommer fc EnBlinl (2004b ) for t he enti re parameter range assuming 



Pfrommer fc Enfilinl ( 12004ai ). Figure 7 for a parametriza- 



tion as adopted in this study. We caution, however, that the minimum-energy condition is 
violated at the outer radio-halo boundary for the range of minimum magnetic-field values 
inferred by this study. 
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candidate is a WIMP. The self-annihilation of WIMPs can produce either monoenergetic 
gamma-ray lines or a continuum of secondary gamma rays that deviates significantly from 
the power-law spectra observed from most conventional astrophysical sources, with a sharp 
cut-off at the WIMP mass. These spectral features together with the expected difference 
in the intensity distribution compared to conventional astrophysical sources allow a clear, 
indirect detection of DM. 



The expected gamma-ray flux due to self-annihilation of WIMPs in a dark-matter halo 
is given by 

» T (An,B)_ <«> <w 



dE 8vrm2 dE 



where (av) is the thermally-averaged product of the total self-annihilation cross section and 



the relative WIM 
per annihilation!^ 



p velocity, m x is the WIMP mass, is the differential gamma-ray yield 
, AQ is the observed solid angle, and J is the so-called astrophysical factor 



- a factor which determines the DM annihilation rate and depends on the DM distribution. 

Given the upper limit on the observed gamma-ray rate, defined as the ratio of the 
event number detected within the observing time T obs , i? 7 (99% CL) = iV 7 (99% CL)/T obs , 
we can place constraints on the WIMP parameter space (m x , (av)). Integrating Eq. ([7j) 
over energy we find 



(av) (99% CL) < i? 7 (99% CL) 



limit , _ 



x 



dE A eS 



dN^E) 



dE 



(8) 



J(Afi) 

where A e g is the effective area of the gamma-ray detector. Because the self-annihilation of 
a WIMP is a two-body process, the astrophysical factor J(AQ) is the line-of-sight integral 
of the DM density squared 



j(ao)= / dn [ dx P 2 x (x,n), (9) 

J AO J 



12 In this work, we have calculated the differential gamma-ray yield per annihilation using 
the Pythia Monte Carlo simulator. 
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where A represents the line of sight. In this work, we have modeled the Coma DM 



distribution with a Navarro, Frenk and White (NFW) profile (INavarro et al 



19971 ). 



Px( r ) = Ps 



-1 



1 + 



-2 



(10) 



where r s is the scale radius and p s is the scale density. Using weak-lensing measurements 



of the virial mass in the Coma cluster and the DM ha 



derive d from iV-body simulation of structure formation (jBullock et al 



o mass-concentration relation 



20011 ) 



Gavazzi et al. 



( 120091 ) find, and list in their Table 1 (note that they define R v i r = R 



100 J) Jw vir 



M< 



200 



9.7(+6.1/ - 3.5) • lO 14 /^ 1 M and C vir = C 2 oo = 3.5(+l.l/ - 0.9), which we translate into 
the density-profile parameters r s = 0.654 Mpc and p s = 4.4 x 10 14 M©/Mpc 3 . Note that 
the uncertainties are not necessarily distributed as a Gaussian, and also arise from the 
choice of dark-matter profile. According to the latest high-resolution DM-only simulations 
of nine rich galaxy clusters, the inn er regions of the smooth density profiles are quite well 



approximated by the NFW formula (IGao et al. 



20121 ) . However, gravitational interactions of 



DM with baryons may modify these predictions. This could give rise to either an increasing 
inner density slope due to adiabatic contraction of the DM component in response to 
cooling baryons in the central regions or a decreasing density slope due to violent baryonic 
feedback processes pushing gas out of the center by, e.g., energy injection through AGNs. 
However, on scales r > 0.45 Mpc or more than 20% of R va (which are of relevance for the 
present work), different assumptions about the inner slope of the smooth DM density profile 
translate to uncertainties in the resulting astrophysical factor. Table [7] lists the astrophysical 
factors calculated for the different VERITAS apertures considered in this work. Table [7| also 
lists the astrophysical factor calculated for the background region, which is used to estimate 
the gamma-ray contamination from DM annihilation in the background region. As long as 
the DM contribution to the event number in the background region is negligible, the upper 
limits derived here directly scale with the astrophysical factor, UL(< av >) oc J -1 . The 
analysis uses a ring region to estimate the background in a ON region. We have to compute 
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the expected level of gamma-ray emission from DM annihilation in the ring region in order 
to check that it is negligible with respect to the level of gamma-ray emission from DM 
annihilation in the ON region. This is equivalent to compute the astrophysical factor of the 
ON and OFF source region since this quantity is related to the rate of DM annihilations. 

The resulting exclusion curves in the ((erf) ,m x ) parameter space are shown in Figure 
[6] for three different DM self-annihilation channels, W + W~, bb, and t + t~. Depending on 
the DM annihilation channel, the limits are on the order of KT 20 to 10" 21 cm s . The 
minimum for each exclusion curve and corresponding DM particle mass is listed in Table 
[HI We stress that these limits are derived with conservative estimates of the astrophysical 
factor J. They do not include any boost to the annihilation rate possibly due to DM 



substructures populating the Coma halo, w hich could sea " 



factor 0(1000) in the most optimistic cases (iPinzke et al. 



e down the present li mits by a 



2011 



Gao et al 



2012 ). 



We also note that when the size of the integration region is increased, the limits on 
(<jv) result from a competition between the gain in the astrophysical factor (J) and the 
integrated background. For integration regions larger than 0.2° in radius, the astrophysical 
factors no longer compensate for the increased number of background events, and the 
signal-to-noise ratio deteriorates. 



6. Discussion and Conclusions 

We have reported on the observations of the Coma cluster of galaxies in VHE gamma 
rays with VERITAS and complementary observations with the Fermi-LAT . VERITAS 
observed the Coma cluster of galaxies for a total of 18.6 hours of high-quality live time 
between March and May in 2008. No significant excess of gamma rays was detected above 
an energy threshold of ~ 220 GeV. The Fermi-LAT has observed the Coma cluster in all-sky 
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survey mode since its launch in June 2008. We have used all data available from launch t o 



2010h . 



April 2012 for an updated analysis compared to published results ( lAckermann et al. 
Again, no significant excess of gamma rays was detected. We have used the VERITAS and 
Fermi-LAT data to calculate flux upper limits at the 99% confidence level for the cluster 
core (considered as both a point-like source and a spatially-extended emission region) and 
for three member galaxies. The flux upper limits obtained were then used to constrain 
properties of the cluster. 

We have employed various approaches to constrain the CR population and magnetic 
field distribution that are complementary in their assumptions and hence well suited to 
assessment of the underlying Bayesian priors in the models. (1) We used a simplified 
"isobaric CR model" that is characterized by a constant CR-to-thermal pressure fraction and 
has a power-law momentum spectrum. While this model is not physically justified a priori, 
it is simple and widely used in the literature and captures som e aspects of more elaborate 



Pinzke &: Pfrommerl (120101 ) 



models such as (2) the simulation-based analytical approach of 
The latter is a "first-principle approach" that predicts the CR distribution spectrally and 
spatially for a given set of assumptions. It is powerful since it only requires the density 
profile as input due to the approximate universality of the CR distribution (when neglecting 
CR diffusion and streaming). Note, however, that inclusion of these CR transport processes 
may be necessary to explain the radio-halo bimodality. (3) Finally, we used a pragmatic 
approach which models the CR and magnetic distributions in order to reproduce the 
observed emission profile of the Coma radio halo. While this approach is also not physically 
justified, it is powerful because it shows what the "correct" model has to achieve and can 
point in the direction of the relevant physics. 

Within this pragmatic approach, we employ two different methods. Firstly, adopting a 
high magnetic field everywhere in the cluster (B 3> -Bcmb) yields the minimum gamma-ray 
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flux in the hadronic model of radio halos which we find to be a factor of 20 (60) below the 
most constraining flux upper limits of Fermi-LAT (VERITAS). Secondly, by matching the 
radio-emission profile (i.e., fixing the radial CR profile for a given magnetic field model) 
and by requiring the pion-decay gamma-ray flux to match the flux upper limits (i.e., fixing 
the normalization of the CR distribution), we obtain lower limits on the magnetic field 
distribution under consideration. Our limits for the central magnetic field range from 
B = 0.5 to 1.4 fiG (for VERITAS flux limits) and from B = 1.4 to 5.5 fiG (for Fermi-LAT 
flux limits). Since all (but one) of these lower limits on Bq are below the values favored 



by Faraday RM, B = 4.7±g;| //G flBonafede et al. 



20101 ) . the hadronic model is a very 



attractive explanation of the Coma radio halo. The Fermi-LAT upper limits start to rule 
out the parameter combination of ag > 0.7 and a > 2.5 for the hadronic model of the 
Coma radio halo. 

Applying our simplified "isobaric CR model" to the most constraining VERITAS 
limits, we can constrain the CR-to-thermal pressure ratio, Xqr, to be below 0.048-0.43 (for 
a CR or gamma-ray spectral index, a, varying between 2.1 and 2.5). We obtain a constraint 
of X C r < 0.1 for a = 2.3, the spectral index predicted by simulations at energies around 220 
GeV. This limit is more constraining by a factor of 1.6 than that of the simulation-based 
model which gives Xcr, < 0.16. This difference is due to the concave form of the simulated 
spectrum which provides more pressure at GeV energies in comparison to a pure power-law 
spectrum of a = 2.3. 

The Fermi-LAT flux limits constrain Xqr to be below 0.012-0.017 (for a varying 
between 2.3 and 2.1), only weakly depending on the assumed CR spectral index. Assuming 
a = 2.1, which is very close to the simulated spectral index for the energy range of 1-3 
GeV, we obtain a constraint which is identical to that from our simulation-based model 
within the virial radius of Xcr < 0.017. That constraint improves to Xqr < 0.012 for an 
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aperture of 0.4° corresponding to a physical scale of R ~ -R200/3 — 660 kpc. These upper 
limits are now starting to constrain the CR physics in self-consistent cosmological cluster 
simulations and cap the maximum CR acceleration efficiency at structure formation shocks 



to be < 50%. Alternatively, this may argue for non-negligible CR transpor ; 



proc esses such 



20121 ). These 



as CR streaming and diffusion into the outer cluster regions (lAleksic et al.l 
are encouraging results in that we constrain the CR pressure (of a phase that is fully 
mixed with the ICM) to be at most a small fraction (< 0.017) of the overall pressure. As 
a result, hydrostatic cluster masses and the total Comptonization parameter due to the 
Sunyaev-Zel'dovich effect suffer at most a very small bias due to CRs. 

We have also used the flux upper limits obtained with VERITAS to constrain the 
thermally-averaged product of the total self-annihilation cross section and the relative 
velocity of DM particles. Modeling the Coma cluster DM halo with a NFW profile we 
derived limits on (av) to be on the order of 10~ 20 to 10~ 21 cm -3 s _1 depending on the 
chosen aperture. These limits are based on conservative estimates of the astrophysical 
factor, where a possible boost to the annihilation rate due to DM substructures in the 
cluster halo has been neglected. Including such a boost could scale down the present limits 
by a factor 0(1000) in the most optimistic cases. 
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Fig. 1. — Left: Smoothed significance map of the Coma cluster calculated from the excess 
VHE gamma-ray events observed with VERITAS over a 4.5° x 4.5° field of view. The color 
scale indicates significance in uni ts of standard deviatio ns. The excess counts were derived 



using a ring-background model ( lAharonian et al. 



2001() . White contours show the X-ray 



counts per second in the 0.1 to 2 keV energy band (8 lev els from 1 to 16 c ts s 1 after 3-pixel 



Gaussian smoothing) from the ROSAT all-sky survey ( jBriel et al. 



19921 ). Right: Same as 



above but with overlaid contours (2 0-180 mJ in 20 mJ steps) from GBT radio observations 



at 1.4 GHz ( jBrown fc Rudnick 



120111 ). where strong point sources have been subtracted. Also 
shown are the 0.2° and 0.4° radii (dashed cyan) considered for the extended-source analyses 
presented here. 
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Fig. 2. — 9 2 distribution from VERITAS observations of the Coma cluster of galaxies. The 
points with error bars represent the ON-source data sample and the filled area is the back- 
ground estimation based on the OFF-source regions. Each bin represents an annulus around 
the Coma cluster core position and the annuli are all of equal area. The data were derived 
from the ring-background model using a 0.2° integration radius. 
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Fig. 3. — Distribution of significances for Figure [D and an integration radius of 0.2°. The 
curve is a Gaussian fit to the data, with mean /i = —0.11 ± 0.0059 and standard deviation 
a = 1.01 ± 0.003, which is consistent with the absence of gamma-ray sources in the field of 
view. 
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Fig. 4. — Integral gamma-ray flux upper limits in narrow energy bands [E ; Ei] for Fermi 
and VERITAS observations of the Coma cluster (given in Tables [2] and E]) for different 
integration radii (arrows with different grey intensities), assuming that a = 2.3 except in 
the 1-3 GeV energy interval where a = 2.1 is adopted. R vu: is the virial radius of the 
Coma cluster, corresponding to 1.25°. These are compared to integrated spectra of the 
same energy interval and aperture, assuming the universal gam ma-ray spectrum of clusters 



(lines with different grey intensities, 



Pinzke fc Pfrommer 



20101 ) . To guide the eye, we show 



the underlying universal integral energy distribution of pion decay gamma-rays, E 1 F 1 {> 
E 7 ), resulting from hadronic interactions of CRs and ICM protons (CR-7r°, dotted). For 
visualization purposes all photon fluxes are weighted with the smallest energy in each interval. 
Note that the Fermi limit for the energy interval of 1-3 GeV within the aperture of 0.4° is 
the most constraining. 
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Fig. 5. — CR-to-thermal pressure ratio, Xcr = (-Pcr) / (-Pth), as a function of radial dis- 
tance from the center of th e Coma cluster. The model predictions, shown by dashed curves 



(IPinzke fc Pfrommer 



2010), have been scaled to match the most constraining VERITAS 



upper limits within 0.2° (solid) and Fermi upper limits within 0.4° (dashed). We compare 
differential X CR profiles (grey) to integrated profiles X CR (< R/ R v i r ) = Pgr dV/ P t h dV 
(black) which we use to compare to the upper limits. 
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Fig. 6. — Limits on the DM annihilation cross section (av) from VERITAS observations 
of the Coma cluster as a function of the DM particle mass m x derived from the VERITAS 
gamma-ray flux upper limits (0.2° aperture) presented in this work. 



-50 - 



Table 1 . Regions of interest in the Coma cluster field of view 



Source RA (J2000) Dec (J2000) 

Core 12 h 59 m 48.7 s +27°58'50.0" 

NGC 4889 13 h 00 m 08.13 s +27°58'37.03' 

NGC 4874 12 h 59 m 35.71 s +27°57'33.37 

NGC 4921 13 h 01 m 26.12 s +27°53'09.59' 



Note. - The cluster core is considered 
both as a point source and a modestly ex- 
tended source. Three central galaxies are 
also considered in point-source searches. The 



choice is based on evidence 



nonth ermal X-ray emission (INeumann et al. 



'or an excess of 



20031 ) at the location of these galaxies. 



Table 2. VERITAS VHE gamma-ray flux upper limits for different regions of interest in the Coma cluster of galaxies 

and surroundings 



Source 


R a [deg] 


N s h 




a 


= 2.1 


Flux U.L. d 
a = 2.3 


a 


= 2.5 


Core 





17 


0.84 


2.59 


(0.78%) 


2.78 


(0.83%) 


2.97 


(0.89%) 




0.2 


-41 


-1.0 


1.96 


(0.59%) 


2.09 


(0.63%) 


2.21 


(0.66%) 




0.4 


-26 


-0.30 


4.44 


(1.3%) 


4.74 


(1.4%) 


5.02 


(1.5%) 


NGC 4889 





3 


0.14 






1.85 


(0.55%) 






NGC 4874 





-14 


-0.71 






1.51 


(0.45%) 






NGC 4921 





-4 


-0.23 






2.41 


(0.72%) 







intrinsic source radius (zero means point source), which is convolved with the 
;amma-ray point-spread function. 



b Net event counts in the source region. 
Statistical significance calculated according to 



Li k Mai (Jl983|) . 



to 



99% confid e nce le vel upper limit in units of 10 8 ph. m 2 s 1 calculated according 



Rolke et al. 



(120051 ) above an energy threshold of 220 GeV, with corresponding fluxes 
in percent of the steady Crab Nebula flux in parenthesis, for different values of the 
spectral index, a. 
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Table 3. Fermi-LAT HE gamma-ray flux upper limits for the Coma cluster core 

1 - 3 GeV 3-10 GeV 10 - 30 GeV 

Spatial Model Flux U.L. a (Significance) Flux U.L. a (Significance) Flux U.L. a (Significance) 

Spectral index a = 2.1 



Point Source 
Disk: r = 0.2° 
Disk: r = 0.4° 



1.882 (0.000) 
2.109 (0.152) 
2.438 (0.201) 



0.759 (0.000) 
0.899 (0.000) 
1.232 (0.619) 



0.671 (0.830) 
0.719 (0.740) 
0.875 (1.387) 



Spectral index a = 2.3 



Point Source 
Disk: r = 0.2° 
Disk: r = 0.4° 



1.946 (0.000) 
2.180 (0.169) 
2.524 (0.246) 



0.788 (0.000) 
0.941 (0.000) 
1.275 (0.742) 



0.667 (0.874) 
0.725 (0.828) 
0.869 (1.390) 



Spectral index a = 2.5 



Point Source 
Disk: r = 0.2° 
Disk: r = 0.4° 



2.008 (0.000) 
2.246 (0.189) 
2.606 (0.291) 



0.816 (0.000) 
0.979 (0.020) 
1.313 (0.856) 



0.663 (0.915) 
0.720 (0.864) 
0.861 (1.387) 
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a 99% confidence level flux upper limit in units of 10 6 ph. m 2 s 



Table 4. Constraints on the CR-to-thermal pressure ratio in the Coma cluster core (simplified, isobaric analytic 
model) for different spatial extensions and predicted fluxes for the energy bands 1-3 GeV and > 220 GeV 

(simulation-based model) 



R a [deg] analytic model: X CR b F^ sim (E) c F UL /F 7iSim (> E) d 
a = 2.1 a = 2.3 a = 2.5 



VERITAS constraints 






0.1 


0.23 


0.97 


1.9 


14.8 


0.2 


0.048 


0.10 


0.43 


2.9 


7.2 


0.4 


0.067 


0.15 


0.62 


4.4 


10.8 



Fermi constraints 






0.035 


0.024 


0.033 


1.4 


1.34 


0.2 


0.024 


0.017 


0.022 


2.1 


1.00 


0.4 


0.017 


0.012 


0.016 


3.2 


0.76 



a Intrinsic source radius (zero means point source), which is convolved 
with the gamma-ray point-spread function. 

b Constraint on the CR-to-thermal pressure ratio, Xqr = 
(Pcr) I (-fth)i which was ass umed to be constant throug hout the cluster 



and calculated according to 



Pfrommer fc Enfflinl (j2004ah . 



by 



Integrated gamma-ray flux from the simulation-based analytic model 



Pinzke fc Pfrommer! (|2010l ): above E = 220 GeV in units of 1(T 9 ph. 

2 



m" 2 s" 1 for VERITAS and for E = 1 - 3 GeV in units of 1(T 6 ph. m 
s _1 for Fermi. 

d Ratio of flux upper limit (-Ftjl) to integrated gamma-ray flux from 
the simulation-based analytic model, with the U.L. based on spectral 
index a = 2.3 in the VERITAS band (E > 220 GeV) and a = 2.1 in the 
Fermi band (E = 1 — 3 GeV) . 



Table 5. Minimum gamma-ray fluxes in the hadronic model of radio halos, where the radio-emitting electrons are 
secondaries from CR interactions, and corresponding minimum CR-to-thermal pressure ratios for Coma 



R a [deg] F 7>min (>£) b 10 4 xX C R, min c 

a = 2.1 a = 2.3 a = 2.5 a = 2.1 a = 2.2, a = 2.5 



VERITAS energy range, E > 220 GeV 






1.6 


0.7 


0.3 


6.7 


6.1 


11 


0.2 


3.1 


1.4 


0.6 


7.8 


7.2 


13 


0.4 


6.3 


2.8 


1.3 


9.8 


9.0 


16 



Fermi energy range, 1-3 GeV 






3.5 


4.8 


6.4 


6.7 


6.1 


11 


0.2 


6.8 


9.3 


12.5 


7.8 


7.2 


13 


0.4 


13.5 


18.6 


25.0 


9.8 


9.0 


16 



intrinsic source radius (zero means point source) , which is convolved 
with the gamma-ray point-spread function. 



b Minimum gamma-ray flux derived from the hadronic model de- 
scribed in § 14.31 Values are in units of 10~ 10 ph. m -2 s _1 for the 
VERITAS energy range and 10 -8 ph. m -2 s -1 for Fermi energy range. 

c Minimum CR-to-thermal pressure ratio, -XcR, m in; in the hadronic 
model described in § 14.31 For simplicity, we duplicate VcR m m for 
the VERITAS and Fermi constraints: for a given realization of the CR 
pressure (and a magnetic field model which is trivial here as we assume 
B S> Bcmb)> we can derive the radio flux as well as gamma-ray fluxes 
in various bands (1-3 GeV, > 220 GeV). 



Table 6. Constraints on magnetic fields in the hadronic model of the Coma radio halo and 
the corresponding CR-to-thermal pressure ratio (at the largest emission radius of 1 Mpc) 
such that the model reproduces the observed radio surface-brightness profile 



Minimum magnetic field, -£>o,min [/wG] a : 
VERITAS constraints Fermi constraints 



a B 


a = 2.1 


a = 2.3 


a = 2.5 


a = 2.1 


a = 2.3 


a = 2.5 


0.3 


0.69 


0.57 


0.48 


1.38 


1.95 


2.68 


0.5 


0.97 


0.80 


0.68 


1.94 


2.74 


3.78 


0.7 


1.40 


1.17 


0.99 


2.80 


3.97 


5.50 






Corresponding 


;X CR (lMpc): 




0.3 


0.46 


1.05 


4.55 


0.11 


0.08 


0.11 


0.5 


0.74 


1.70 


7.47 


0.18 


0.13 


0.17 


0.7 


1.09 


2.59 


11.55 


0.27 


0.19 


0.26 



a The parameters of the magnetic field are the magnetic de- 
cline, ob, and the central field strength, B , which are defined 
by B(r) = B [n e (r)/n e (0)] as . In all cases, we used the most 
constraining R = 0.4°; see § 14.41 for details. 



Table 7. Astrophysical factors 



R [deg] (J) signal [GeV 2 cm- 5 sr] a (J) bkg [GeV 2 cm- 5 sr] 



5.7 x 10 16 1.3 x 10 14 (negligible) 

0.2 8.1 x 10 16 4.4 x 10 14 (< 0.01 (J) signal , negligible) 

0.4 9.4 x 10 16 1.3 x 10 15 (~ 0.01 (J) si „ nal , negligible) 



Note. — (J)bkg * s ^ ne astrophysical factor calculated for the background 
region (ring method) and is used to estimate the level of gamma-ray con- 
tamination from DM annihilation, a is the size ratio of the ON and OFF 
source regions. 
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Table 8. Upper limits on the DM annihilation cross section times velocity (av) from 
VERITAS observations of the Coma cluster. 



Channel 


R [deg] 


m x [GeV] 


(av) [cm 3 s x ] 


W+W" 





2000 


1.1 x 10~ 20 




0.2 


1900 


4.3 x 10" 21 




0.4 


1900 


8.4 x 10~ 21 


bb 





3500 


1.2 x 10~ 20 




0.2 


3400 


4.4 x 10~ 21 




0.4 


3500 


8.7 x 10~ 21 


r + T - 





670 


2.4 x 10~ 21 




0.2 


650 


9.1 x 10" 22 




0.4 


660 


1.8 x 10~ 21 



Note. - Upper limits are shown for different 
integration regions and DM particle mass m x , and 
are derived from the VERITAS gamma-ray flux 
upper limits presented in this work. 



